Correlation Effects and the High-Frequency Spin Susceptibility of an Electron Liquid: 

Exact Limits 
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Spin correlations in an interacting electron liquid are studied in the high-frequency limit and 
in both two and three dimensions. The third-moment sum rule is evaluated and used to derive 
exact limiting forms (at both long- and short-wavelengths) for the spin-antisymmetric local-field 
factor, lim^^oo G- (q, uj). In two dimensions limo,-,^ G_ (q, uj) is found to diverge as 1/q at long 
wavelengths, and the spin-antisymmetric exchange-correlation kernel of time-dependent spin density 
functional theory diverges as 1/q 2 in both two and three dimensions. These signal a failure of the 
local-density approximation, one that can be redressed by alternative approaches. 
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I. INTRODUCTION 



The detailed incorporation of exchange and correlation 
(xc) into the dynamic response of an interacting electron 
liquid remains a challenge to theory. The mean field ap- 
proach of the random phase approximation (RPA), which 
neglects xc effects, becomes increasingly inadequate as 
both the dimensionality and density of electrons is low- 
ered. Sum rules provide exact constraints on the response 
of an electron liquid and thus are extremely useful as 
tests of any scheme that purports to include xc effects. 
Alternatively they can be used to constrain particular 
approximations such as those embodied in local field fac- 
tors G± (q, uj) . The static behavior of both the spin- 
symmetric (+) and spin-antisymmetric (— ) local field 
factors has been well studied^ and various fitting formu- 
lae have been proposed that conform to the known exact 
constraints 2,3 - 4 . Their dynamic forms have attracted a 
great deal of recent attention since they exhibit far richer 
behavior than their static counterparts. 



However, up until now, most dynamical studies have 
focussed on the spin-symmetric case; in this paper we re- 
port on a derivation of the exact behavior of the dynamic 
spin-antisymmetric response in the high-o) limit within 
the linear response scheme. In particular, the first and 
third-moment sum rules are derived (Seem and then 
used in Sec lIIII to determine the exact asymptotic be- 
haviour of lim^^oo G- (q, uj) in the low- and high-g limit. 
The xc kernel in time-dependent density functional the- 
ory is closely related to these local field factors and in 
Sec lIVI it will be shown that the high-frequency spin- 
antisymmetric xc kernel must then diverge in the long 
wavelength limit. This divergence demonstrates that the 
dynamical xc potential in time-dependent spin-density 
functional theory (TDSDFT), as applied to a homoge- 
neous system, must be a nonlocal functional of the local 
density. 



II. MOMENT SUM RULES 

To describe the formalism needed to derive the mo- 
ment sum rules of an electron liquid we will first need 
the electron particle-hole operator pk a (q), defined as 



Pko-(q) = 4 CT c k+q!(T , 



(1) 



where c qo . and c qCT are the fermionic creation and anni- 
hilation operators, respectively, and a denotes the spin 
projection. The standard one-particle density operators 
are then 



/5<r(q) = ^k<r(q), P(q) = ^ /Wq): ( 2 ) 

k 

and, for the spin equivalent, 



krr 



-1 for <7 =j, 
-1 for <t =| . 



(3) 

In the thermodynamic limit the Hamiltonian for a d- 
dimensional system of N electrons immersed in a uniform 
neutralizing background of volume £1 is given by (fi = 1), 

H = 5> (q)c qCT c qa + ^ E nq)[^P(q)p(-q) - 1] 



qtr 



q/0 



where Wo(q) = q 2 /2m, and V(q) is the Fourier trans- 
formed Coulomb potential, i.e. V^q) = 4ire 2 /q 2 for 
d = 3, and V(q) = 2ire 2 /q for d — 2. From linear re- 
sponse theory the retarded charge-density \c and spin- 
density xs responses are given by 



xc(q,*) 



[p(q,*),|S(-q,O)])0(t), 



(5) 



and 



Xs(<l,t) = -([§(<l,t),s(-q,0)])9(t), (6) 



where 0(t) is the Heaviside step function and () denotes 
averaging over the equilibrium ensemble as specified by 



2 



the Hamiltonian in Eq. (@J. Note that to obtain a phys- 
ical spin-spin response from xs(q, t) in Eq. © a factor 
g 2 p B /4 is required where g is the gyromagnetic ratio and 
Hb is the Bohr magneton. 

As usual, the spectral representation of xio.^), the 
Fourier transform of x(q> can be obtained from 



x(q,w) 



eL/ Imx(q, w') 

7T a/ — U3 — M + 



(7) 



Analyticity of xsio., w ) m the upper half of the complex- 
ly plane requires that Imxs (q, be an odd function of 
u) and so, in the the high-w limit, xg(q, expands as 



lim Xs(q,w) 



1=1 



(21-1) 



,21 



(8) 



where the moments are given by 

( W (0) = ^ ^^Im Xs (q, W ), 

s(q, i),s(-q, o) 



5 



9t 



(9) 



The first moment yields the well known /-sum rule, 

= 2n^ (q) (10) 

where n = AT/ fi is the macroscopic homogeneous particle 
density. The derivation of the /-sum rule invokes number 
conservation via the continuity equation, 



dt 



= -iq-j<r(q,*) 



(ii) 



where j<x(q, t) is the current density fluctuation operator 
for particles with spin a 



j<r(q,£) = — y^( k + ^q)pk<r(q^)- 



(12) 



Thus number density is separately conserved for electrons 
with spin up and down since there is no term in the equi- 
librium Hamiltonian Eq. Q which can flip spin. That 
correlation effects play no role in the first moment may 
be understood by noting that in the high-w limit elec- 
trons which interact via a velocity-independent potential 



cannot be influenced by the effects of others when the 
time interval is infinitesimal. 

The third-moment sum rule^, after straightforward but 
lengthy calculations, gives&£ 

(lu^) = 2n Wo (q)[ W 2(q) + u Q (q) (12/ d)(E KE ) +7 d (q)] 

where (Eke) is the average kinetic energy per electron 
in the interacting system and id(q) is given by 

J*(q) =-£ V(k)(q • k) 2 [5(q - k) - S(k)], (14) 
where for q ^ 



5(q) - AT^qM-q)), 
5(q) - iV- 1 {p(q) / 5(-q)) ! 



(15) 
(16) 



define the static magnetic and usual structure factors, 
respectively. For completeness, we express the structure 
factors in terms of the spin-resolved pair correlation func- 
tions, 



S(k)-1 = - I ^r( m (r)- Pn (r))expHk-r), (17) 



and 



S(k) - l =zJ dr GmW + 9n(r) - 2) exp(-tk • r). 

,( 18 ) 

where g-\-\(r) and gn(r) are the spin-parallel and antipar- 
allel contributions, respectively, to the electron pair cor- 
relation function g(r), i.e. 



9(r) = q[ffn( r ) +-9T1 Ml- 



(19) 



It has been noted^ that the third-moment sum rule dif- 
fers for the charge and spin density response, in contrast 
to the first moment sum rule which is obeyed by both \c 
and xs- The physical root of this qualitative difference 
can be inferred from the continuity equation for longi- 
tudinal spin-current, used in the derivation of Eq. \YA\ . 
namely 



*^q-ja(q,«) = X>°(q + k) - Wo (k)] 2 pk CT (q,t) + -J2v(p)^-cl a (t)p(p,t)c^ p+cia (t). (20) 

k kp 

I 

Hence, it can be shown that the presence of the inter- across spin-up and spin-down electrons. Thus, although 
electron potential allows a transfer of momentum density the total momentum density must be conserved, the mo- 
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mentum density of each spin species is noA This is why 
the third-moment sum rule, essentially an expression of 
momentum and number density conservation, must differ 
for charge and spin density susceptibilities. 

An expression for I c i (q) has been known for some time 
for three-dimensional systems^, and here we just record 
the limiting expressions, 



and 



/s(q-0) 



7 3 (q -> oo) = 



^(q) 



<4(q) 



[i-m(o)], 



[1-2^(0)], 



(21) 



(22) 



where w p (q) = [2nwo(q)^(q)] 1/ ' 2 is the plasma frequency. 
Note that in Eqs. (|21I - I22|) we have applied^ # TT (°) = °> 
and that the high-q limit is to be taken such that <C 
w(q) <C lo. 

In two dimensions we find 

/ 2 (q) = ^ / dkk{R(c l ,k)[S(k)-l]-k[S(k)-l]j 

(23) 

with 



i?(q,k) 



dO q 2 + k 2 cos 2 6 + 2kq cos 9 



io n (q 2 + k 2 + 2kqcos9) 1 / 2 
The limiting behavior is then 



(24) 



ia(q-»0) = — / dkk 2 [S(k) - S(k)}, (25) 



2m 



I 2 (q^oo) = --Wp(q)# n (0). 



III. LOCAL-FIELD FACTORS 



(26) 



Another way to express the correlation effects in the 
spin response is to introduce the notion of a local-field 



factor, G_(q, w), which modifies the effective field felt 
by a single electron, as expressed by 



xs(q^) 



xo(q> w ) 



l + F(q)G_(q,cj)x (q,cj)' 



(27) 



Here x (q, cj) is the modified Lindhard function for the 
interacting system, i.e. where exact interacting occupa- 
tion numbers, n(k), are used, namely 



xo(q> w ) = o 



n(k) - n(k + q) 



-y- 

Q ^ uj + w (k) - tu (k + q) + i0+ 



(28) 



The n(k) therefore replace the usual non-interacting 
fermion occupation number, n(k), in the original Lind- 
hard function, xo(q> w). Since the exact occupation num- 
bers are generally unknown we follow Ref. |10| and intro- 
duce a separate local factor, G„ to write the modified 
Lindhard function in terms of the original Lindhard func- 
tion XO: i- e - 



Xo(q> w ) 



xo(q,w) 



1 + F(q)G n (q,u;)xo(q,w) 



(29) 



Thus we have 



xs(q,w) = - 



xo(q,w) 



l + V(q)G_(q,w)xo(q,w) 



(30) 



where G_ (q, w) = G_ (q, w) + G n (q, w). 

The third-moment sum rule can again be utilised to 
determine the high-a; limit of the local-field factors, i.e. 
lim^^oo G(q, u>) — G^°°'(q). First we will require the 
high-cj expansion of the Lindhard function, i.e. 



lim Xo(q,w) 



2nwo(q) 2ncj 2 (q) 



12 

^o(q) + -t(E K e)o 
d 



(31) 



where (Eke) a is the average kinetic energy per electron in a non- interacting system. After insertion of Eq. (|31|l into 
the high-w expansion of the spin-density response Eq. (|30() we observe that manifestations of the local field factor 
appear at 0(l/w 4 ), i.e. 



lim xs(q,w) 



2nwo(q) 2no>Q(q) 



wc(q) + ^-(Eke)o - y(q)2nG (o ° ) (q) 
d 



O 



(32) 



r 



By comparison with the first- and third-moment sum rules, Eq. JTUl and Eq. JT^, we infer that limiting form 
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of the local-field factor must be given by 



G^q) 



^(q) 



G 



dnV(q) 



((Eke)o-(Eke))- (33) 



The kinetic part of the above expression gives a diverging 
contribution in the large q limit which, previous authors 8 
have argued, must be absorbed by a proper treatment 
of the self-energy of the proper polarizability function, 
i.e. the modified Lindhard function. The introduction of 
the additional local- field factor, G n , allows us to demon- 
strate this explicitly, and indeed a high-a; expansion of 
Eq. H28JI together with the definition of the local field- 
factor, Eq. (J20J), gives 



^(q) 



<> 



dnV(q) 
leaving the desired result, 



G^q) 



((Eke)o - (Eke)), (34) 



Mq) 



(35) 



This result can also be derived from the work of 
Zhu and Overhauserii based on the equation-of-motion 
approach 1 —. In three dimensions the limiting forms are 



G^Cq^O) 
G^q-oo) 



[25(0) - 1], 
[4ff(0) - 1], 



(36) 
(37) 



Note that if the negative sign of the local-field factor at 
high to, by Eq. (|36fl . persists at finite ui then we have 
the possibility that a pole in xs(q, ^) exists 1 — for real 
ui and sufficiently large r a . Physically, this collective 
mode corresponds to long- wavelength spin- wave modes, 
undamped and oscillating at frequencies given by the so- 
lution to 



l + y(q)G_(q,o;)xo(q,w) = 



(38) 



In two dimensions the limiting forms of the high- 
frequency local-field factors are 



G^q^O) 



G^q-oo) = 5(0). 



1 



Airnq 



dkk 2 [S(k) - S(fc)],(39) 
(40) 



IV. DISCUSSION 

The physical concept behind the introduction of the 
local-field factor is identical in spirit to the methodology 
of DFT in which xc effects are incorporated in an effective 
local field experienced by each electron. In fact, within 
the frequency-dependent linear response formalism 14 , the 
unpolarized spin-antisymmetric local-field factor is sim- 
ply related to the xc energy functional i£ xc [n] via, 



G_(q,w) 



KU (q, ui) - Kll (q, ui) _ K~ (q, ui) 



2V(q) 



2V(q) 



(41) 



where K™ (q, ui) is the Fourier-transform of the xc ker- 
nel of TDSDFT. 



S 2 E xc [n] 



Sn l7 (r,t)Sn iyl (r,t') 



(42) 



From Eq. 1361) and Eq. H39|) we now note that 
lim q ^o lim^^oo K~ c (q, ui) diverges as q~ 2 in both two 
and three dimensions in notable contradistinction to 
K£ c (q, ui) which is well-behaved in the same limit. This 
striking behavior of the spin-antisymmetric exchange- 
correlation kernel does not by itself lead to any concern 
since the corresponding macroscopic quantities (i.e. xs 
as given by Eq. (|30l) ) remain well-defined in the low-q 
and high-w limit and it follows that the observable (i.e. 
induced magnetization) also remains well-defined. Phys- 
ically, the divergence in momentum-space implies that 
K~ c (r, ui), and hence the xc hole, is of infinite extent, 
i.e. it decays very slowly, especially so in two dimen- 
sions. The divergence is then a signature of a transition 
from a bounded xc kernel to an unbounded one. Thus 
a strict local-density approximation (LDA) 14 for the xc 
potential becomes untenable in time-dependent systems, 
at least in the high-frequency regime. In other words, 
the dynamical local field cannot be specified by the local 
density alone and other approaches to TDSDFT are re- 
quired. This divergence is particularly potent given that 
it occurs even in a homogeneous system, as emphasised 
by Vignalei^ who very recently reported similar conclu- 
sions about the dynamic spin response at finite-a->. 

One possible remedy is via the weighted density 
approximation 1 ^ which retains nonlocality and has been 
demonstrated to give more accurate results than func- 
tional approximations involving local-density contribu- 
tions. Indeed it has been arguediZ*^ that such an ap- 
proach should be able to recover at least part of the 
the full Van der Waals interaction, the origin of which is 
clearly due to dynamical correlation between electrons, 
and where it is well known that the conventional LDA 
fails. Another approach which overcomes the limitations 
of the LDA but nonetheless still permits a local descrip- 
tion is current density functional theory 19,20 where the 
local current density is required in addition to the local 
density to describe time-dependent inhomogeneous sys- 
tems. 

Finally, we note that many-body calculations based on 
Hedin's GW approximation^ take correlations into ac- 
count directly without invoking the local mean-field ap- 
proximations of DFT, and have recentl}^ been demon- 
strated to yield accurate ground-state energies of inho- 
mogeneous systems where LDA fails drastically. How- 
ever, such many-body approaches are numerically inten- 
sive and it remains to be seen whether they can be suc- 
cessfully applied to excited states within feasible compu- 
tational effort. 
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